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We carried out light curve solutions of the ultrashort-period binaries with MS components observed by Kepler. All six 
targets turned out almost in thermal contact with contact or slightly overcontact configurations. Two of them, KID 4921906 
and KID 6309193, are not eclipsing but reveal ellipsoidal and spot variability. One of the components of KID 8108785 
exhibits inherent, quasi-sinusoidal, small-amplitude variability. KID 12055255 turned out a very rare case of ultrashort- 
period overcontact binary consisting of two M dwarfs. Our modeling indicated that the variability of KID 9532219 is due 
to eclipses but not to S Set pulsations as it was previously supposed. 
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1 Introduction 

Most of the contact binaries consisting of solar-type com¬ 
ponents have orbital periods within 0.25 d < P < 0.7 d. 
Rucinski (1992) found that they show a short-period limit 
at about 0.22 d. But the modern large surveys allowed to 
discover targets beyond this limit (Maceroni & Rucinski 
1997; Maceroni & Montalban 2004; Weldrake 2004; Rucin¬ 
ski 2006, 2007; Pribulla et al. 2009; Lohr et al. 2014; Qian 
et al. 2014, etc.). 

The well-studied binaries with periods around and be¬ 
low cut-off limit (ultrashort-period systems) and MS com¬ 
ponents are quite rare but important objects for the astro¬ 
physics, especially for the stellar evolution. 

Until two decades ago the statistics of contact binaries 
was particularly affected by the trend to discover relatively 
large light-curve amplitudes, but very recently, the huge sur¬ 
veys (ROTSE, MACHO, ASAS, SuperWASP, etc.) led to 
discoveries of many low-amplitude systems. 

Lately, Kepler provides unique, unprecedented precise, 
extended and nearly uninterrupted data set for a large num¬ 
ber and variety of stars. Above two thousands eclipsing bi¬ 
naries (EBs) were identified and included in the Kepler EB 
catalog (Prsa et al. 2011, further PRSA11; Slawson et al. 
2011). According to the initial classification the sample con¬ 
sists of 1261 detached, 152 semidetached, 469 overcontact 
binaries, 137 ellipsoidal variables (they are not EBs!), and 
146 uncertain systems. The automated fitting of the light 
curves by a polynomial chain and the artificial intelligence 
based tier EBAI (Prsa et al. 2008) were used for the morpho- 
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logical classification and determination of the global param¬ 
eters of some systems of the Kepler EB catalog. 

Several individual Kepler EBs were studied in details 
(e.g., Southworth et al. 2011; Steffen et al. 2011; Welsh et 
al. 2011; Winn et al. 2011). This paper is devoted to standard 
(’’manual”) light curve solutions of the ultrashort-period bi¬ 
naries from the Kepler EB catalog. 

2 Selection of targets and preliminary 
preparation of the data 

As a continuation of our interest in binaries with extremely 
small orbital periods (Dimitrov & Kjurkchieva 2010), we 
reviewed the Kepler EB catalog. It contains 22 targets with 
periods P < 0.23 d (for better statistics we put an upper 
limit 0.23 d that is slightly bigger than the cut-off limit of 
0.22 d). 

Although Prsa et al. (2011) reported about visual in¬ 
spection of every target and culling about 27% of their ini¬ 
tial sample as nonEBs (RR Lyr-, 7 Dor-, S -Set type stars), 
we carried out new visual review of the chosen 22 targets. 
We concluded that 13 ultrashort-period stars are S Set vari¬ 
ables due to the shape of their curves (KID 10417986, KID 
8912468, KID 8758716, KID 10855535, KID 9612468, 
KID 8912730, KID 9898401, KID 7375612, KID 5872696, 
KID 10684673, KID 6699679, KID 6287172, KID 
11825204). Our classification of these stars is supported by 
their high temperatures (above 6000 K). 

One ultrashort-period star, KID 9472174, turned out to 
be the known sdB+dM binary 2M1938+4603. 

The shapes of the light curves of two ultrashort-period 
targets, KID 8288741 and KID 12350008, seem as those of 
detached eclipsing systems (with almost flat maxima) but 
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their amplitudes of 1-2 mmag and periods are absolutely 
inappropriate for such configurations. Our preliminary at¬ 
tempts to reproduce their data required enormous big third 
light. This means that KID 8288741 and KID 12350008 are 
not the true variable stars. It is meaningless to model these 
data because the obtained parameters would not be reliable. 
Matuevic et al. (2012) noted that the automated classifica¬ 
tion and modeling of Kepler EBs, even with locally lin¬ 
ear embedding, failed for small amplitudes (^0.001 mag). 
KID 8288741 and KID 12350008 are such cases but the rea¬ 
son to be removed from the further analysis is the discrep¬ 
ancy between the light curve shape appropriate for detached 
eclipsing binary and the negligible amplitude of variability 
(smaller even than those of ellipsoidal variations). 

As a result of the selection procedure we gathered sam¬ 
ple of six ultrashort-period binaries: KID 4921906, KID 
1572353, KID 8108785, KID 6309193, KID 12055255, 
KID 9532219. More conclusive confirmation of their du¬ 
plicity would be spectroscopic data. The objects would re¬ 
quire 8m-class telescope due to low brightness and short pe¬ 
riod at the same time. Another possibility for proper classifi¬ 
cation is the color-period diagram (Duerbeck, 1997) but the 
presence of third component or close neighbor can alter the 
color information. So we will use the important criterion for 
sure classification of the targets: the successful reproducing 
of their light curves by effects of duplicity. 

Table 1 presents information for our 6 targets from the 
Kepler EB catalog: Kepler magnitude epoch To of 
the primary minimum; orbital period P\ mean temperature 
T m ; metallicity [Fe/H]\ radius R\ log g\ reddening E(B-V); 
type of configuration; ratio T 2 /Ti of the temperatures of the 
components; fillout factor FF ; crowding (a number between 
0 and 1 that specifies the fraction of flux in the optimal pho¬ 
tometric aperture due to the target star with respect to the 
total flux from all sources). 

More than 50000 points are available for each target in 
the Kepler archive. The detailed review of these data re¬ 
vealed that the shape of the light curves almost did not 
change during the different observational quarters. However 
the levels of maxima and minima and the amplitudes of the 
light curves of the most targets changed considerably during 
the different quarters causing the big thickness of the total 
light curves (from all quarters). This effect is probably due 
to problems of the satellite steering and automated reduc¬ 
tion of the data. That is why we assumed that it is reasonable 
to model data from several consecutive quarters rather than 
all data simultaneously. Moreover, we established that it is 
more appropriate to use the raw Kepler data and to make 
own reduction and de-trending (Dimitrov et al. 2012). 

We chose to model the data from quarters Q0 and Q1 
for all targets excepting KID 9532219 which observations 
began later. Their folded light curves are presented in Figs. 
1 - 6 . 
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Fig. 1 The light curve of KID 4921906 and our fit (top) 
and the corresponding residuals (bottom) 
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Fig. 2 The light curve of KID 1572353 and our fit (top) 
and the corresponding residuals (bottom) 

3 Light curve solutions 

The shapes of the light curves (Figs. 1-6) implied that our 
targets are nearly contact or overcontact systems that was 
expected for their short orbital periods. The depths of the 
two minima of the most curves are close that means that the 
binaries are almost in thermal contact. The small light am¬ 
plitudes of some targets point at low orbital inclinations and 
consequently, considerable contributions of the ellipsoidal 
variations. 

Taking into account the foregoing considerations we 
carried out modeling of the photometric data by the code 
PHOEBE (Prsa & Zwitter 2005) using the following proce¬ 
dure. 
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Table 1 Parameters of the targets from the EB catalog 


Kepler ID 

m K 

[mag] 

T 0 -2400000 
[days] 

P 

[days] 

Tm 

[K] 

[Fe/H] 

R 

[*G>] 

log# 

E(B-V) 

Type 

t 2 /t 1 

FF 

crowd 

4921906 

15.203 

54965.148704 

0.213732 

4800 

-0.394 

0.99 

4.41 

0.109 

OC 

0.958 

1.033 

0.648 

1572353 

15.234 

55001.306290 

0.228900 

4470 

-0.261 

0.77 

4.54 

0.076 

OC 

0.883 

1.051 

0.9 

8108785 

14.757 

54964.641700 

0.228840 

4354 

-0.298 

0.67 

4.61 

0.053 

OC 

0.855 

0.414 

0.921 

6309193 

13.674 

55002.434080 

0.212500 

5367 

-0.383 

1.09 

4.38 

0.098 

OC 

1.012 

1.01 

0.92 

12055255 

15.866 

54964.966610 

0.220950 

4093 

0.106 

0.73 

4.50 

0.066 

9 

? 

? 

0.566 

9532219 

16.118 

55001.945246 

0.198155 

5031 

-0.539 

1.06 

4.38 

0.159 

9 

? 

? 

? 
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Fig. 3 The light curve of KID 8108785 and our fit (top) 
and the corresponding residuals (bottom) 


Initially the primary temperature was fixed as Ti=T m . 
The mean target temperatures T m (Table 1) have been es¬ 
timated using dedicated pre-launch ground-based optical 
multi-color photometry plus Two Micron All Sky Survey 
(2MASS) J, K , and H magnitudes (Skrutskie et al., 2006). 

We varied the secondary temperature T 2 , mass ratio q, 
orbital inclination i and potentials S 2 i 5 2 (and simultaneously 
relative radii 7 * 1,2 and fillout factors FFi^)- In order to re¬ 
produce the O’Connell effect we added cool spots on the 
primary and varied their parameters: longitude /3, latitude 
A, angular size a and temperature factor k = T sv jT\. 

We adopted coefficients of gravity brightening 0.32 and 
reflection effect 0.5 appropriate for late stars. The limb- 
darkening coefficients were chosen according to the tables 
of Van Hamme (1993). 

Finally, we varied the component temperatures around 
T m to search for the best fit. 

The derived parameters of the targets corresponding to 
our light curve solutions are given in Table 2 while Table 
3 presents information about the spot parameters. The for¬ 
mal PHOEBE errors of the fitted parameters are quite small 
(Tables 2-3) due to the high precision of the Kepler data. 

It should be noted that the ephemerides from the Kepler 
EB catalogue (PRSA11) did not lead to the best fits. For 
this aim we varied the parameter ’’shift” of PHOEBE which 
final values are given in column 2 of Table 2. 
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Fig. 4 The light curve of KID 6309193 and our fit (top) 
and the corresponding residuals (bottom) 



Fig. 5 The light curve of KID 12055255 and our fit (top) 
and the corresponding residuals (bottom) 
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Table 2 Parameters of the best light curve solutions 


KID 

shift 

i 

[°] 

9 

,, 

T 2 /T 1 

t?l,2 

FF 1j2 

„mean 
r l,2 

h 

config 

variability 

4921906 

0.0126 

21.04±0.02 

0.5977±0.0005 

4528±14 

0.960 

3.0634±0.0002 

0.0 

0.424 

0.66 

C 

ellipsoidal 





43483=12 


3.06343=0.0002 

0.0 

0.334 






(26,4) 

(0,47) 


(0.958) 





(OC) 


1572353 

-0.01923 

43.84±0.02 

0.9577±0.0003 

45253=20 

0.945 

3.6883=0.0003 

0.0 

0.383 

0.58 

c 

grazing 





42773=13 


3.688±0.001 

0.0 

0.374 



eclipses 



(33) 

(0.29) 


(0.883) 





(OC) 


8108785 

0.00147 

61.78±0.02 

0.83093=0.0003 

47113=4 

0.901 

3.46643=0.0004 

0.0 

0.395 

0.62 

c 

eclipses 





4246±3 


3.46643=0.0006 

0.0 

0.360 






(69.9) 

(0.33) 


(0.855) 





(OC) 


6309193 

-0.0196 

13.445±0.014 

0.8007±0.0003 

54353=2.5 

0.874 

3.4165 ±0.0005 

0.0 

0.398 

0.68 

c 

ellipsoidal 





4751 ±2 


3.4165 ±0.0005 

0.0 

0.359 






(21) 

(0.60) 


(1.012) 





(OC) 


12055255 

0.0128 

39.79±0.05 

0.96373=0.0005 

3701±8 

1.011 

3.6326±0.0022 

0.128 

0.392 

0.55 

OC 

grazing 





3744±10 


3.6936±0.0012 

0.019 

0.377 



eclipses 



(?) 

(?) 


(?) 





(?) 


9532219 

0.00647 

68.42±0.06 

0.9320±0.0022 

4950±22 

0.983 

3.5687±0.0022 

0.172 

0.398 

0.55 

OC 

eclipses 





48673=21 


3.5873±0.0028 

0.164 

0.388 






(?) 

(?) 


(?) 





(?) 



Note: C means contact binaries, OC - overcontact systems 
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Fig. 6 The light curve of KID 9532219 and our fit (top) 
and the corresponding residuals (bottom) 


Table 3 Parameters of the spots and third light of the tar¬ 
gets 


KID 

P 

[°] 

A 

[°] 

a 

[°] 

k 

h 

4921906 

75.0±0.1 

88.1±0.1 

11.0±0.1 

0.901±0.001 

0 

1572353 

91.5±0.1 

329.9±0.1 

17.6±0.1 

0.900±0.002 

0 

8108785 

126.3±0.2 

279.6±0.3 

24.9±0.2 

0.939±0.001 

0 

6309193 

63.9±0.1 

338.4±0.1 

10.3±0.2 

0.940±0.001 

0 

12055255 

59.6±0.2 

89.9±0.2 

14.9±0.2 

0.900±0.001 

0 

9532219 

81.5±0.6 

70.09±0.04 

10.1±0.7 

0.914±0.007 

0.78 


The synthetic curves corresponding to the parameters 
of our light curve solutions are shown in Figs. 1-6 as red 
continuous lines. 

The derived photometric mass ratio requires a special at¬ 
tention. Usually this parameter needs spectral observations. 
But the spectral mass ratio for short-period binaries with fast 
rotation is quite imprecise due to the line blending while the 
geometrical effects and mass ratio are much more important 


in defining the shapes of their light curves and the eclipse 
depths. Hence, the values of the photometric mass ratio for 
ultrashort-period Kepler binaries should be considered with 
high confidence. 

There are differences between our solutions and those 
obtained by automated fitting cited in the Kepler EB catalog 
(the values in parentheses in Table 2). 

(1) The most orbital inclinations in PRSA11 are bigger than 
ours. 

(2) The mass ratios in PRSA11 are considerably smaller 
than our values (Table 2). 

(3) We had to assume a third light / 3 only for KID 9532219 
(Table 3) while the values of this parameter in PRSA11 
were nonzeros for the all six targets. We found a weaker 
star close to KID 9532219 (at distance around 6 arcsec) 
that could explain the presence of / 3 in the light curve 
solution of this target. 

(4) There are no values of the potentials Di ? 2 (and thus of 
the relative radii 7 * 1,2 and fillout factors FFi^) as well 
as of relative luminosities in PRSA11 that means that 
the automated modelling leads to approximate results. 

There are several explanations for the differences be¬ 
tween our ’’manual” solutions and those of the automated 
fitting. 

(a) The EBAI method fixes / 3 = 1 - crowding while this 
parameter is invoked in our procedure only when a sat¬ 
isfactory fit is impossible for any combination of the pa¬ 
rameters. The bigger value of Z 3 requires bigger light 
amplitude and correspondingly bigger orbital inclina¬ 
tion. Moreover, due to the correlation between third 
light and mass ratio for partially eclipsing binaries the 
bigger third light leads to the larger mass ratio for fixed 
inclination angle. 

(b) Although all light curves exhibit O’Connell effect (Figs. 
1-6) the EBAI method does not take into account the ef¬ 
fect from spots. Then the results of the automated mod¬ 
elling should be assumed as first approximations of the 
fitted parameters. 
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4 Analysis of the results 

The analysis of our light curve solutions of the ultrashort- 
period binaries from the Kepler data base led to several con¬ 
clusions. 

(1) The configurations of all targets are almost the 
same: 4 of them are contact binaries (KID 4921906, 
KID 1572353, KID 8108785, KID 6309193) and 2 
are slightly overcontact systems (KID 12055255 and 
KID 9532219). 

(2) Two targets, KID 4921906 and KID 6309193, are 
noneclipsing systems. The registration of their low- 
amplitude ellipsoidal and rotational (due spots) vari¬ 
ability is possible due to the high Kepler accuracy. 
There are quite many objects of this type in the Ke¬ 
pler EB catalog and hence, it should be called catalog 
of photometric ally-double stars. 

(3) To model the data of KID 9532219 we had to assume 
considerable contribution of a third light Z 3 . It originates 
probably from the close neighbor. We do not exclude 
this fainter star to be the true variable because the size of 
the Kepler pixels projected onto the sky (4 arcsec) cou¬ 
pled with the high star density near the Galactic plane 
lead to a non-negligible likelihood of associating an EB 
event with the wrong star (PRSA11). 

(4) The target KID 12055255 is very rare case of ultrashort- 
period binary consisting of two M dwarfs. The system 
seems to be overcontact. We are not aware of any over¬ 
contact binary consisting of M dwarfs. 

(5) Feldmeier et al. (2011) classified the target 
KID 9532219 as a S Set variable (from their ground- 
based photometric observations of the Kepler FoV for 
the pre-launch survey). The Kepler data and our mod¬ 
eling certainly reveal that this target (or its neighbor) is 
eclipsing binary. 

( 6 ) Recently KID 6309193 was noted as ’’false” in the Ke¬ 
pler EB catalog. Our modeling showed that it is not 
eclipsing star but its light curve can be reproduced by 
ellipsoidal variations of binary. However we do not ex¬ 
clude the alternative the observed variability to be inher¬ 
ent to some of its two considerably fainter neighbors (at 
distances around 4.5 arcsec and 7 arcsec). 

(7) The ripples in the residuals of KID 8108785 (Fig. 3) 
seem to be quasi-sinusoidal with period of 55 min and 
amplitude of 0.005 mag. Although this variability is vis¬ 
ible during all phases, it could be explained by inherent 
changes of some of the components of this system due 
to its small orbital inclination. 

Although all considered light curves reveal O’Connell 
effect and are reproduced by cool spots we have not stud¬ 
ied long-term light curve changes as a result of the evolu¬ 
tion of these surface inhomogeneities. To do this one should 
have detailed information about the problems of the satellite 
steering which introduce different trends into the data. 


5 Conclusions 

Our investigation added six new members to the small fam¬ 
ily of the studied ultrashort-period binaries with MS com¬ 
ponents. 

All six targets turned out to be contact or slightly over¬ 
contact configurations which components are almost in ther¬ 
mal contact. Two targets, KID 4921906 and KID 6309193, 
are not eclipsing but reveal ellipsoidal and spot variabil¬ 
ity. KID 12055255 is very rare case of ultrashort-period 
binary consisting of two M dwarfs. We established that 
KID 9532219 is an eclipsing binary but not a S Set vari¬ 
able as it was previously supposed. One of the components 
of KID 8108785 exhibits inherent, quasi-sinusoidal, small- 
amplitude variability. 

The analysis of the results showed that reliable global 
parameters of the Kepler binaries might be obtained only 
by standard method of light curve solutions. The automated 
EBAI method models the Kepler data quite roughly. Obvi¬ 
ously, the numerous and exclusive precise Kepler data de¬ 
serve precise light curve solutions by the standard way. 
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